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Abstract. We present a high-resolution 5-GHz radio con- 
tinuum map of the starburst galaxy NGC 2146 made with 
MERLIN and the VLA (A-array), in a search of radio su- 
pernovae and supernova remnants expected to be already 
produced by the most massive stars in the starburst. At 
5 GHz, about 20 point sources were detected earlier by 
Glendenning & Kronberg (1986) in the central 800 pc of 
NGC 2146. Our observations with higher sensitivity and 
resolution made with MERLIN and the VLA confirms the 
detection of 18 sources, and resolves 7 of them. Additional 
1.6-GIIz MERLIN observations disclose 9 sources coinci- 
dent in position with those detected at 5 GHz, which al- 
lows us to derive their spectral indices af g, (S^ ~ i^"). 
Only 3 sources have indices (a < 0) consistent with syn- 
chrotron emission from supernova remnants or radio su- 
pernovae, while the others have very steep inverted spectra 
(a > 0). 

We suggest that the sources with positive spectral in- 
dex are optically thick ultra-compact and/or ultra-dense 
H II regions with high electron densities and high emission 
measures (EM > 10^ cm~^ pc). Minimum energy require- 
ments indicate that these regions may contain up to 1000 
equivalent stars of type 06, comparable to the number 
of stars found in super starclusters. When compared with 
M82, the galaxy NGC 2146 lacks however a large num- 
ber of supernova remnants. We suggest that NGC 2146 is 
experiencing a burst of star formation stronger than that 
in M82, but being in a younger phase. In this phase, only 
few stars have already exploded, whereas the others cause 
strong thermal emission from compact, optically thick ion- 
ized gas regions, around the young super starclusters. 
We may, however, not exclude an alternative scenario in 
which strong free- free absorption at 1.6 GHz in foreground 



ionized gas with very high emission measures (EM > 10* 
cm^^ pc) hides a certain number of supernova remnants, 
thus rendering for some sources the observed inverted 
spectra. 

Key words: Hil regions - ISM: supernova remnants - 
Galaxies: individual: NGC 2146 - Galaxies: starburst - 
Radio continuum: ISM 



1. Introduction 

NGC 2146 is a starburst galaxy located at a distance of 
14.5 Mpc (l"w 70 pc). Its infrared luminosity measured 
with IRAS at 60 ^m and 100 ^iul is 6.6 x 10^° Lq. This 
value and the large 25 ^m to 60 ^m flux ratio place the ob- 
ject in the lower part of the IR luminosity/spectral index 
plane populated by superluminous galaxies without active 
nuclei (Hutchings et al. |1990|). NGC 2146 has a central 



molecular ring (Jackson & Ho 1988, Young et al. 1988), 
and an outflow of hot gas along the minor axis driven by 
supernova explosions and stellar winds in the starburst 



region (Armus et al. 1995, Delia Ceca et al. 1999). These 
characteristics reveal a strong similarity to the prototype 
starburst galaxy M 82 (at a distance of 3.2 Mpc and with 
an infrared luminosity of 2.4 x 10^° L0), although NGC 
2146 does not have a companion that may have triggered 



the starburst (Fisher & TuUy 197(:). It has been suggested 
that the starburst in NGC 2146 is the result of a far 



1988 



evolved merger (Condon et al. 1982; Young et al 
Hutchings et al. 1990D , but a fully convincing kinematic 



and material trace of the merger has not yet been found. 
The nuclear region is partly obscured by a strongly absorb- 
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ing dust lane (Benvenuti et al. 1975 ). The high extinction 
(Av > 5 ma g; You ng et al. 1988 ; Hutchings et al 



1990 



Smith et al. 1995| ) makes the optical observation of the 
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inner star-forming regions impossible. 
Strong non-thermal radio emission from the centre of 
NGC 2146 has been detected by Kronberg & Bicrmann 
( |1981| ; hereafter KB) by C ondon et al. ( |1982D , and by 
Lisenfeld et al. ( [1996| , |1997D . About 20 point sources were 
detected at 5 GHz in the central 800 pc of NGC 2146 



by Glendenning & Kronberg (1986) using the VLA, and 
these were interpreted as radio supernovae (RSN) or su- 
pernova remnants (SNRs). In this paper we report the 
detailed spatial distribution of these point sources, to- 
gether with their number and fluxes derived from radio 
continuum observations at 5 GHz obtained with MER- 
LIN and the VLA. From additional 1.6 GHz MERLIN 
observations we obtained the spectral index of 9 of these 
sources, which allow to identify their nature. Our results 
open the possibility of a direct comparison with similar 
ob servat ions of M82 (Kronberg et al. 



1985; Muxlow et 



al. 1994 ) and other strong and nearby starburst galaxies, 
such as NGC 1808, NGC 4736 and NGC 5253 (Saikia et al. 



199C; Duric & Dittmar 1988; Turner et al. 1998) 



2. The observations and image processing 

NGC 2146 was observed for 23 h with MERLIN (6 an- 
tennas) in November 1996 (one day) and February 1997 
(two days). The observing frequency was 4.994 GHz (A 
= 6 cm), with a bandwidth of 15 MHz in both circular 
polarizations; the data were taken in spectral-line mode 
(15 X 1-MHz channels). OQ208 (2.39 Jy) and 0602+780 
(0.126 Jy) were used as flux and phase calibration sources, 
respectively. 

NGC 2146 was observed in the L-Band with MERLIN (6 
antennas, including the Lovell telescope) in May 1997, for 
a period of 15 h. The observing frequency was 1.658 GHz 
(A = 18 cm), with a bandwidth of 15 MHz in both circu- 
lar polarizations. The passband and relative gains of the 
antennas were determined using the point source calibra- 
tor 0552-1-398, with a flux density of 2.31 Jy; the compact 
source 0602-1-780 (0.109 Jy) was used to determine the 
telescope phases. 

Moreover, we used observations of NGC 2146 taken in 
September 1985 with the VLA |^ (27 antennas) in the A, 
B and C array configurations; the observing frequency was 
4.885 GHz, with an integration time of 18 h. 
The 5 GHz data from MERLIN and the VLA A-array 
were combined to increase the brightness sensitivity and 
uv coverage, in order to improve the information about the 
weaker sources already marginally detected in the 5 GHz 
MERLIN observation. To combine the two data sets it was 
necessary to shift the MERLIN phase reference position 
slightly. A re-weighting was also necessary to compensate 
for different origins of the data sets. Several images were 
produced using the AIPS task IMAGR, and deconvolved 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 



with the CLEAN algorithm (Hogbom |1974D . The details 
of these maps such as restoring beam, the applied weight- 
ing, and the rms noise, are summarized in Table The 
rms noise in source-free areas of the images is consistent 
with the expected thermal noise levels. 
Since the Wardle-telescope was not present in the obser- 
vations, the length of the shortest baseline of the MERLIN 
array at 1.6 GHz is 9.2 km (Darnhall-Lovell telescopes). 
At this frequency the interferometer is then not sensi- 
tive to structures larger than ■^l^l'^^^ = 4" (280 pc). At 
5 GHz the shortest baseline of the MERLIN array is 9.5 
km (Darnhall-Mark2 telescopes), so that -d^^J^^ = 1.3" 
(90 pc). The combination of the 5 GHz MERLIN data 
with the VLA data gives a shortest baseline of 0.68 km, 
so that i?^Sx' = 18" (1250 pc). The short-spacing effect 
does not influence the results at either frequency since the 
size of the sources we are searching for is smaller than 
1" (100 pc). Furthermore, most of the sources are either 
unresolved or have extensions of a few beam widths and 
hence the flux measurement is not influenced by the lack 
of short spacings. 

3. The results 

Figure |l| shows the uniformly weighted ABC-array VLA 
map (A) of NGC 2146. Due to the higher resolution (0'.'37 
X 0'.'32 instead of 1'.'2 x 1'.'2) the comparison of this map 
with that of KB reveals a larger number of compact radio 
sources 18 instead of 5). Any feature above 0.1 mJy 
(3cr level) in the 5 GHz MERLIN + VLA map (C) is taken 
to be a source. Our map lacks the long 'tongue' of emis- 
sion to the SE, seen in the map of KB. The global features 
found by KB, such as the S-shape of the central emission 
region (by KB assumed to be a 'triple' source, but now re- 
solved into a larger number of sources), its position angle 

o o 

(~ 105 ), and the main 'bar' of radio emission (~ 128 ), 
are consistent with those in our map. Convolution of our 
VLA map to the 1'.'2 resolution of the KB map shows 
that the peak brightness is of the order of 6 mJy/beam 
as measured by KB, and that the SE tongue of emission 
is still present, though less extended than in the earlier 
maps. From the map shown in Fig. ^ it is evident that 
the point sources present in the SE part of the galaxy are 
more numerous than in the NW part (13 compared to 4, 
assuming that the centre of symmetry of NGC 2146 is near 
the strong central source 37.6-1-24.2). The extended emis- 
sion of the SE region is also more pronounced, giving the 
impression of an asymmetry in NGC 2146, in favour of a 
stronger starburst activity along the SE arm. ^^CO inter- 
ferometer observations (Young et al. 1988) show a higher 
concentration of molecular gas in this part of the galactic 
disk. 

Table |2] contains the positions and the fluxes of the eigh- 
teen compact sources derived from our observation, listed 
by coordinate name with increasing R .A., f ollowing the 
convention of Kronberg & Wilkinson (1975). The refer- 
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Map 
code 


Frequency 
(GHz) 


Observation 


Weighting 


Beam 
HPBW(") 


Noise 
(pijy beam~^) 


A 


5 


VLA 


Purely uniform* 


0.37 X 0.32 


~ 20 


B 


5 


MERLIN 


Purely uniform* 


0.04 X 0.04 


~ 50 


C 


5 


MERLIN + 


Uniform 


0.15 X 0.14 


~ 35 






VLA A-array 








D 


1.6 


MERLIN 


Purely natural** 


0.19 X 0.15 


~ 35 



' : using the robust parameter 



using the robust parameter = +5 [AIPS task IMAGR]. 



ence position is at a = OO'^IS^OWO and S = 78°21'00'.'0 
(J2000). The coordinates of the emission peaks were de- 
termined from a Gaussian fit of each component using 
the AIPS task JMFIT (which also gives the absolute po- 
sitional uncertainties) in the uniformly weighted 5 GHz 
MERLIN + VLA map (C). From this map we also derived 
the peak and integrated flux densities at 5 GHz (S5). The 
same quantities were derived at 1.6 GHz (Si.g) from the 
naturally weighted MERLIN map (D). The numbers in 
parentheses are the errors of the derived quantities. Only 
nine sources are detected at 1.6 GHz, due to a lack of sen- 
sitivity in the MERLIN observations, and possibly also to 
the steep positive spectral indices which may place the re- 
maining sources below the 3 a level. The spectral indices 
(af g) were obtained from the integrated flux densities at 
both frequencies, where we adopt the convention . 
The flux densities have been computed using the AIPS 
task JMFIT, when the sources were either unresolved or 
shell-like, and using the task BLSUM for sources with 
extended structures impossible to fit by a single Gaus- 
sian. The integrated fiux densities were calculated up to a 
3(7 level to avoid confusion with the extended background 
emission. 

In Table ^ we report the deduced quantities for the eigh- 
teen sources: the total emitted power at 5 GHz (L5) and 
1.6 GHz (Li e), for a distance of 14.5 Mpc, based on the 
peak flux values of Table ||, and the brightness temper- 
ature calculated from the integrated flux intensities at 
5 GHz (Table ^). The diameter (D) was meas ured as 
described in Appendix 1 of Burns et al. ( 1979| ), where 

R.„ 



R 



Tcsolvcd 



= 2 



(0.85 • HWHMbo 



Runrcsoived = 0.85 • HWHMbcam, Rmap the radius taken 
from the map, and HWHMbcam the half width at half 
maximum of the Gaussian beam. The diameter for the 
unresolved sources is 8.6 pc (Table ||) and their brightness 
temperatures are obviously lower limits because their true 
dimensions are not known. 



3.1. Sources with negative spectral indices 

3.1.1. The "central source" 37.6+24.2 

37.6+24.2 has a large negative spectral index of g = 
-1.1, consistent with optically thin synchrotron emission. 
The source has a position offset of ~ 1'.'15 north-east of the 
dynamic center of NGC 2146 (R.A. = 06^18™3?;4 ± 0!4, 
Dec. = 78°21'23'.'2 ± 2") deduced from the 1^00(1 - 0) 
rotation curve measured with the Plateau de Bure inter- 
ferometer, adopting a systemic velocity of 850 km s'^ . The 
source is not resolved in the 0'.'15 resolution 5-GHz map 
(C), and its brightness temperature is too high to be due 
to free-free emission from ionized gas. The strength of the 
source allowed us to map it with a resolution of 0'.'04 us- 
ing only the MERLIN observations at 5 GHz (B). Even in 
this case the source is not resolved, indicating that it is a 
very strong and compact object (less than 3 pc diameter). 
This source was considered earlier to be the true nucleus 
of NGC 2146 (KB). However, the properties of the source 
are consistent with those found for RSN and SNRs, and 
because of its steep spectrum and sma ll size , we believe 
that this source is a RSN (Weiler et al. 1986 ). 



3.1.2. Source 38.9+22.5 

38.9+22.5 has also a large negative spectral index af g = 
-1.1. The object is resolved both at 5 GHz and 1.6 GHz, 
and has the structure of a partially filled shell of ^ 30 
pc diameter (Fig. |^). As a consequence of its morphology, 
the size obtained from a Gaussian fit corresponds to its 
largest angular extent, which may explain the relatively 
low brightness temperature at 5 GHz. Similar to the pre- 
vious source we believe that also in this case we are dealing 
with either a SNR or a RSN. While the steep spectrum 
indicates it to be a RSN, its large size and the shell-like 
morphology agree better with that of a SNR. 
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Table 2. Observed quantities of the compact sources 



R.A. Dec. Peak Flux Integrated Flux Peak Flux Integrated Flux Spectral 
(J2000) (J2000) Ss Ss Si.e Si.e Index 

Source Name 06''18'" 78°2l' (mjy/beam) (mjy) (mjy/beam) (mjy) ale Notes^^ 



35.0+28.0 


34! 96 


28'.' 


0.27 


(0.04) 


0.27 


(0.06) 










u 


35.3+30.9 


35^32 


30'.' 9 


0.32 


(0.03) 


0.75 


(0.11) 










r 


36.6+27.5 


36! 62 


2/.' 5 


0.56 


(0.04) 


0.61 


(0.06) 


0.15 


(0.04) 


0.15 (0.07) 


+1.3 


u 


36.8+24.0 


36! 78 


24'.' 


0.50 


(0.03) 


2.39 


(0.20) 










r 


37.6+24.2 


37!59 


24'.' 2 


1.56 


(0.04) 


1.56 


(0.06) 


3.90 


(0.04) 


5.30 (0.08) 


-1.1 


u 


37.7+23.7 


37!69 


23'.' 7 


0.31 


(0.03) 


0.45 


(0.08) 


0.26 


(0.04) 


0.26 (0.07) 


+0.5 


u 


37.9+23.7 


37!91 


23'.' 7 


0.35 


(0.06) 


1.42 


(0.09) 


0.30 


(0.04) 


0.30 (0.07) 


+1.3 


r 


38.1+26.1 


38! 12 


26'.' 1 


0.15 


(0.04) 


0.15 


(0.06) 










u 


38.2+24.3 


38! 21 


24'.' 3 


0.22 


(0.03) 


0.29 


(0.07) 










u 


38.4+23.8 


38! 44 


23'.' 8 


0.20 


(0.03) 


1.36 


(0.25) 










r 


38.6+24.0 


38! 59 


24'.' 


0.28 


(0.03) 


0.61 


(0.11) 










r 


38.9+22.5 


38! 90 


22'.' 5 


0.30 


(0.03) 


1.15 


(0.20) 


0.29 


(0.03) 


4.10 (0.46) 


-1.1 


r* 


39.0+17.9 


39! 05 


17'.' 9 


0.28 


(0.04) 


0.28 


(0.06) 


0.16 


(0.04) 


0.16 (0.07) 


+0.5 


u 


39.2+20.7 


39! 16 


20'.' 7 


0.26 


(0.04) 


0.26 


(0.06) 










u 


39.5+17.7 


39! 48 


ir.'7 


0.28 


(0.03) 


0.71 


(0.11) 


0.27 


(0.04) 


0.36 (0.08) 


+0.6 


r 


40.2+18.1 


40! 16 


18'.' 1 


0.23 


(0.03) 


0.37 


(0.08) 


0.18 


(0.04) 


0.18 (0.07) 


+0.6 


u 


41.2+17.3 


41! 16 


17'.' 3 


0.17 


(0.04) 


0.17 


(0.06) 










u 


41.4+15.0 


41!36 


15'.'0 


0.25 


(0.04) 


0.25 


(0.06) 


0.50 


(0.04) 


0.50 (0.07) 


-0.6 


u 



r: resolved only at 5 GHz; r*: resolved at both frequencies; u: unresolved 
Table 3. Diameters, luminosities, and brightness temperatures of the observed sources. 



Ls Brightness Li.e 

D from Peak Flux Temperature^ from Peak Flux 

Source Name (pc) x 10^^ erg • s"^ • Hz'^ (K) x 10^^ erg • s • Uz'^ 



35.0+28.0 


8.6 


6.78 


>1260 




35.3+30.9 


9.4 


8.04 


3010 




36.6+27.5 


8.6 


14.07 


>2840 


3.77 


36.8+24.0 


18.4 


12.56 


2455 




37.6+24.2 


8.6 


39.20 


>7265 


98.00 


37.7+23.7 


8.6 


7.79 


2095 


6.53 


37.9+23.7 


24.8 


8.80 


820 


7.54 


38.1+26.1 


8.6 


3.77 


> 700 




38.2+24.3 


8.6 


5.53 


> 1350 




38.4+23.8 


22.6 


5.03 


945 




38.6+24.0 


9.1 


7.04 


2610 




38.9+22.5 


28.5 


7.54 


500 


7.29 


39.0+17.9 


8.6 


7.04 


>1305 


4.02 


39.2+20.7 


8.6 


6.53 


> 1210 




39.5+17.7 


9.4 


7.04 


2850 


6.78 


40.2+18.1 


8.6 


5.78 


>1725 


4.52 


41.2+17.3 


8.6 


4.27 


> 790 




41.4+15.0 


8.6 


6.28 


>1165 


12.56 



^ Rayleigh- Jeans 
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36.6+27.5 












Fig. 1. The 5 GHz pure uniformly weighted A,B,C-array VLA map (A), with a beam (HPBW) of 0'.'37 x 0'.'32. The 
peak flux density is 1.6 mJy/beam and the contour interval is 0.08 mJy/beam (4 a), up to 40% of the peak surface 
brightness. The smaller maps are sub-images of 9 sources from the MERLIN + VLA A-array uniformly weighted map 
(C, resolution 0'.'15 x 0'.'14). For these sources we have calculated spectral indices. In 37.6+24.2 the contour interval 
is 0.1 mJy/beam (3 a) with the first contour at 0.15 mJy/beam (4 a). In 36.6+27.5, 39.0+17.9, 39.5+17.7, 40.2+18.1, 
and 41.4+15.0 the contour interval is 0.1 mJy/beam (3(t) with the first contour at 0.1 mJy/bcam (3 a). In 37.7+23.7, 
37.9+23.7, and 38.9+22.5 the interval is 0.05 mJy/beam (2 a) with the first contour at 0.15 mJy/beam (4 a). 



3.1.3. Source 41.4+15.0 



3.2. Sources with positive spectral indices 



This source has a negative spectral index of g 



0.6, 



again consistent with synchrotron emission. It is not re- 
solved and the derived brightness temperature of 1200 K is 
a lower limit. The spectral index suggests that this source 
is a SNR. 



There are 6 sources which have a positive spectral index 
(af 6 > 0). Four of them (36.6+27.5, 39.0+17.9, 40.2+18.1 
and 37.7+23.7) are not resolved, so that their brightness 
temperatures are only lower limits. The other two sources 
(37.9+23.7 and 39.5+17.7) are resolved into amorphous 
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Fig. 2. Left panel: source 38.9+22.5 at 5 GHz; the contour interval is 0.05 mJy/beam (2 a) with the first contour at 
0.15 mJy/beam (4 a). Right panel: source 38.9+22.5 at 1.6 GHz; the contour interval is 0.05 mJy/beam (2 a) with 
the first contour at 0.1 mJy/beam (3 a). 



structures. The brightness temperatures, though relatively 
high, are consistent with free-free emission. 

4. Discussion 

In the following Sects. (4.1 and 4.2) we try to give an 
explanation of the physical nature of the thermal and non- 
thermal sources detected in our observations. 

Thermal sources (positive spectral index) 

A positive spectral index with a slope higher than 0.5 
is expected to arise from free-free emission either in the 
turnover region, where r « 1, or in the optically thick 
region, where r > 1. The free-free optical depth in the 
radio domain is 



r = 8.2 X 10" 



GHz 



-2.1 



EM 



pccm 



-1.35 



(1) 



(Mezger & Hend erson |l967t Carlstrom & Kronberg |1991 
Wills et al. 1998 ), where v is the observing frequency, Tc 
the electron temperature, and EM the emission measure 
defined as 



EM = 



(2) 



where Ue is the electron density and L the pathlength. Us- 
ing in Eq.(l) the electron temperature Tc = 8000 K de- 
rived from radio recombination line observations of NGC 



2146 (Zhao et al. 199q), we obtain for r = 1 the emis- 
sion measure YM].li « 6 x 10'' cm^^ pc at 1.6 GHz 
and EM^ GHz « 7 x 10^ cm"^ pc at 5 GHz. For sources 
with a very steep positive spectral index these values is 
only a lower limit since for these sources r > 1. If the 
thermal gas is optically thin (r w 0.1) and has an electron 
temperature of 8000 K, the applied source detection limit 



at 5 GHz of 0.1 mJy corresponds for an unresolved source 
to an EM of ~ 4 x 10'^. Using the value EM(t=1), and for 
the pathlength L the diameter D of the sources (Table ^ , 
we derive the electron densities given in Table [J. 
If we assume that the free-free emission is thermal 
bremsstrahlung from Hii regions photoionized by hot, lu- 
minous stars, it is possible to estimate the Lyman con- 
tinuum photon rate Nly [s~^] from the measured radio 
flux, and from NLy it is possible to estimate the number 
of exciting stars. For optically thin emission, Hjellming & 
Newell ( |1983D give 



n: 



thir 

Ly 



= 3.1 X lO"* 



f To 



V4000K 



-0.16 



V4.9GHzy 



ISOpcy V7mJy 



(3) 



where d is the distance of th e galaxy; for optically thick 
emission, Turner et al. (1998) give 



N^l'JJ"'' = 47r • 



R 



Vcm" 



(4) 



where q;b is the coefficient for recombinations to all but 
the ground state of hydrogen (=3.2 x 10""'^'^; tabulated for 
different electron temperatures and densities in Hummer 



& Storey 1987). Having thus computed the Lyman con- 
tinuum photon rate NL y, we use the compilation of stellar 
properties by Panagia ( 1973] ) to derive the number of stars 
embedded in the individual objects observed by us. The 
number of equivalent 06 stars^ is given in Table |^. 
Are there objects in NGC 2146 which have emission mea- 
sures of the order of EM w 7 x 10^ cm^^ pc, and if so, can 



This is an arbitrary, though often used, unit. Using instead 
tlie also common unit of 07V or 05.5 stars will change the 
number of equivalent stars in Table Wlby ± 50 %, respectively. 
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there be stellar clusters of 1000 to 2000 06 stars associ- 
ated with them? 

'Classical' Galactic Hii regions have EMs smaller than 10^ 
cm~^ pc. Applying the definition of the emission measure, 



Carlstrom & Kronberg (1991) derived a hypothetical up- 
per limit of 4.5 x 10^ cnT^^pc for the EMs of the ionized 
gas in M 82, assuming that the gas is uniformly distributed 
throughout a 500 pc disk, with an average density less 
than 300 cm^-^ (Duffy et al. |1987 ). F rom radio and mil- 



limeter o bserva tions (Seaq uist e t al 
Kronberg [l99ll Wills et al. |l99j) 



1985 



Carlstrom & 
values of up to 5 x 10^ 
cm"'' pc have been derived for M82. 
The EMs required to explain the observations of NGC 
2146 are higher and suggest that the sources we are deal- 
ing with must have outstanding properties. Large EMs 
of up to 10^ cm~^ pc have been found in our galaxy in 
compact and ultracompact Hii regions (see Churchwell 
1990| ) with typical sizes of 0.01 pc. High EMs (10^-10^ 



cm pc) have been recently also found in two galaxies. 



the starburst dwarf galaxy NGC 5253 (Turner et al. |1998D 
and the blue compact galaxy Henize 2-10 (Kobulnicky & 
Johnson 2000 ) , in less compact though ultra-dense H ii re- 
gions with size between 0.1 and 8 pc. The parameters of 
our unresolved thermal sources (turnover frequency, size, 
EM), together with their estimated equivalent number of 
06 stars, are consistent with these sources. As far as the 
resolved thermal sources (37.9-1-23.7; 39.5-1-17.7) are con- 
cerned, the equivalent number of 06 stars seems clearly 
unrealistic, especially for 37.9-1-23.7, which is surely seen 
in the optically thick regime ^ The star cluster associated 
with 37.9-1-23.7 is not seen at optical wavelengths, nor is 
any other super star cluster. 

We thus propose that the objects with positive spectral 
index are in reality a mixture of thermal and non-thermal 
sources, too compact to be fully resolved in our map. How- 
ever, at the single frequency of 5 GHz it is difficult to sepa- 
rate between thermal and non-thermal emission. Further- 
more free-free absorption by ionized gas could play a role 
in obscuring non-thermal sources, such as SNRs and/or 
RSN, embedded in or located behind optically thick Hii 
regions (with different sizes and densities). Of course such 
an explanation presents difficulties in being physically re- 
alistic. When we assume that the turnover frequency is 
equal to or larger than 5 GHz, the EM of the thermal 
sources should be larger than 10* cm~^ pc, implying H ii 
regions that surround the non-thermal objects, with very 
small sizes (0.01-0.1 pc), but containing hundreds of mas- 
sive stars. If the thermal sources are instead located in 
front of, and are not associated with the SNRs or RSN, 
then the number of stars should be even larger. However, 
the statistical probability of a line-of-sight coincidence of 



^ The sources with less steep spectra belong to the turnover 
region, so the equivalent number of 06 stars should lie between 
the values computed for the optically thick and optically thin 
case (which is a much smaller number). 



SNR/RSN and non- associated compact H il regions should 

be very small for the six observed sources. 

Fortunately, the discovery of very dense thermal sources 



by Kobulnicky & Johnson (2000), with sizes up to 8 pc and 
containing ~ 750 07V equivalent stars, makes the possi- 
bility less extreme of actually dealing with thermal sources 
of parsec scales with a very high stellar density. In addition 
to this, we cannot exclude that the free-free absorption is 
affecting only the 1.6 GHz observations. In this case the re- 
quired EM would be ~ lO"^ cm~^ pc, and as a consequence 
the sizes of the H ii regions could be larger, and the stellar 
density smaller. Better constraints on the exact turnover 
frequency, the mixture of thermal and non-thermal emis- 
sion, the Lyman continuum rate from the free-free emis- 
sion and the equivalent number of 06 stars present both, 
in the resolved and unresoved sources, can be obtained 
from high-resolution observations at still higher frequen- 
cies (for instance 8 and 15 GHz). At these frequencies 
the free-free radio emission will originate either from opti- 
cally thin H II regions (in this case the spectral index will 
change to a slope of ~ -0.1), or from non-thermal sources 
no longer affected by free-free absorption (in this case the 
slope will be typical for synchrotron emission). 
The presence of a number of compact and ultra-compact 
Hii regions in NGC 2146 has also been claimed by Puxley 



et al. (1991) and Zhao et al. (1996) using observations of 
H53a (43 GHz) and H92a (8 GHz) recombination lines. To 
explain both, the strength of the 8 GHz continuum and the 
recombination lines, they require a two-component model 
of Hii regions with high electron densities 10^ cm^^) 
and sizes of 0.2 pc, and lower electron density H ii regions 



(- 10^ 



cm ■^) with parsec-scale sizes. 



The reason why we do not see such ultra-compact/ultra- 
dense sources in the starburst galaxy M 82, which in many 
respects is similar to NGC 2146, is unknown. A sequence 
for the starburst development has been proposed for some 
starburst galaxies, including M82, by Rieke et al. (1988) 
and this picture may also apply to NGC 2146. Accord- 
ing to this sequence, NGC 5253 for instance, which has a 
flat radio spectrum and similar compact sources as NGC 
2146, is classified to be in a very young phase of a star- 
burst, while M82 is showing an older phase in which the 
massive stars have already evolved off the main sequence, 
with many SNRs clearly visible. NGC 2146 is probably 
in between these two phases because it still has optically 
thick Hii regions like those in NGC 5253, and a smaller 
number of SNRs than M82, but nevertheless, more simi- 
lar to M82, exhibits a global radio spectral index with a 
non-thermal slope of -0.74. 

4-2. Non-Thermal sources 

4.2.1. Luminosities and energetics 

Earlier ra dio fre quency observations of NGC 2146 (Mc- 
Cutcheon 1973 KB) show that the galaxy is filled 
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Table 4. Sources with positive spectral index 



Source 
Name 



Electron Density 
„™-3 



Equivalent number 
of 06 Stars^ 



Equivalent number 
of 06 stars^ 



36.6+27.5 
37.7+23.7 
37.9+23.7 
39.0+17.9 
39.5+17.7 
40.2+18.1 



> 2800 

> 2800 

> 1600 

> 2800 

> 2700 

> 2800 



2060 
2060 
16750 
2060 
2410 
2060 



1310 
960 
3040 
600 
1520 
790 



^ For optically thick emission (Eq. 
^ For optically thin emission (Eq. t 



with non-thermal synchrotron emission of spectral index 
a w —0.74^ This situation resembles that of M82, where 
the spectral index is a « —0.6 (White & Becker 1992; 
Kiihr et al. |1981 ). Following current theories, this syn- 
chrotron emission is due to SN, RSN, and SNRs which 
preferentially occur in the starburst region, as supported 
by the observed X-ray emission and outflows of hot gas 
along the minor axis. Using minimum-energy arguments 
(Miley |1980D , Muxlow et al. ( |1994D have computed the 
number of SNRs 2000) necessary to explain the ob- 
served flux density at 5 GHz of the extended radio com- 
ponent in the central part of M 82 and consequently the 
time (4 x 10^ yr) necessary to produce such a flux den- 
sity, assuming a supernova rate of 0.05 yr~^. Following 
t his c omputation, we use the relation given by Longair 
(|1994D to calculate the minimum energy requirements for 
synchrotron radiation: 



8 X 106.(1 + /3)4/7. 

ycm-'/ 



VHz/ 



erg • s 



Hz- 



4/7 



[ergj 



(5) 
(6) 



where (3 is the energy ratio of relativistic protons to elec- 
trons, taken to be 1, V is the volume of the considered 
region, v is the frequency and h^, the luminosity. The 5- 
GHz flux of the central region of NGC 2146 (assumed to 
be a disk with a radius of 750 pc and thickness of 100 
pc) is 0.25 Jy, which corresponds to a minimum energy 
of ~ 5 X 10^'^ erg. To produce this energy, ~ 10^ evolved 
SNRs similar to those in the LMC (~ 5 x 10^^ erg each) 
are required. A supernova rate of 0.15 yr~^ (see Sect. 



* This value has been confirmed using the total flux den- 
sity of 475 ± 39 mjy obtained from new single-dish 5 GHz 
observations of NGC 2146, obtained with the 100m Effelsberg 
radio telescope of the MPIfR, and the total flux at 20 cm of 
1228 mJy taken from th e Gre en Bank 1.4 GHz Northern Sky 
Survey (White & Becker |l99j). 



4.2.2) could produce the extended non-thermal radio emis- 
sion over a time of 6 x 10** yr. The individual events are 
probably too old and faint still to be seen individually, 
but their number and/or strength must be higher than 
in M82, making the small number of SNRs detected sur- 
prising. The strength of the overall magnetic field is not 
taken into account, though it must play an important role 
in increasing/decreasing the synchrotron emissivity and 
the global spectral index of a galaxy. 
We also computed the minimum energies and the mag- 
netic field strengths of the three SNR candidates presented 
in the previous section. The equation for the magnetic field 
(Longair 1994 ) is 



B„i„ = 9.3 X 10-^ 



V 



2/7 



.1/7 



[Gauss] 



(7) 



with the same units as in Eq. ^. The results are compiled 
in Table p| . Compared with the values obtained by Duric et 
al. (1988[) for non-thermal compact sources in NGC 4736, 



the agreement of minimum energy requirements and the 
higher magnetic field strength again favour an identifica- 
tion with fairly young SNRs. The total emitted power at 
5 GHz (Table ^ is typically more than ten times higher 
than those found by Duric et al. ( 1988| ) and than that of 
Cas A. Instead it falls in the upper end of the range of the 
M 82 sources identified as young SNRs or RSN. 



4.2.2. The number of SNRs and RSN in NGC 2146. 

The star formation rate (SFR) of NGC 2146, calculated 
using the extinction-corrected Ha fiux, is ~ 6 Mq yr~iand 



is three times higher than that of M 82 (Young et al. 1988). 
The same result is obtained from the FIR fiux (Moshir et 



al. 1990), using the relation derived by Thronson & Tele- 



sco (1986) and Thronson et al. (1991). The values for the 
SFRs in NGC 2146 and M82 then are 15 Mq yr-^_and 
5 M0yr^i respectively. Van Buren & Greenhouse ( 1994 ) 
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Table 5. Sources with negative spectral index 



5. Conclusions 



Source 


Uniin 


Bniin 


Name 


(10^" erg) 


(10~^ Gauss) 


37.6+24.2 


1.5 


410 


38.9+22.5 


6.0 


135 


41.4+15.0 


0.4 


215 



use the FIR flux to calculate the supernova rate in M 82 
to be 0.06 yr"^ Applying this method to NGC 2146 we 
obtain a supernova rate of 0.15 yr~^, which is again higher 
than the one in M82. From these results, and despite the 
fact that the volume of the starburst region in NGC 2146 
is larger than that in M 82, we can infer that the bursts in 
both galaxies have different strengths. As a consequence, 
we should have a higher number of SNRs and RSN in NGC 
2146 than in M82. 

However, when we place the ~ 40 SNRs and RSN de- 
tected inM 82 at 5 GHz (Kronberg et al. |l985l Muxlow 



et al. |1994[ ) at the distance of NGC 2146, given the sen- 
sitivity (~ 0.1 mJy at 3cr level) of our MERLIN + VLA 
map (C) we expect to find at least 8 sources . The num- 
ber of sources detected in our observations is only one 
third of that expected. For 9 sources, however, we do not 
have spectral indices, and their continuum spectra may 
be distorted by free-free or synchrotron self-absorption as 
explained above. A larger number of SNRs than those de- 
tected in NGC 2146 would also be expected if, as argued 
by Armus et al. ( |1995[ ), SNRs are the possible cause of 
the hard X-ray component in NGC 2146 We may there- 
fore speculate that the starburst in NGC 2146 is young, 
or repetitive, so that we are presently observing a strong 
and young burst still lacking a large number of SNRs and 
RSN, but instead emitting strong optically thick thermal 
emission from ionized gas regions, possibly associated with 
super starclusters. A first, also strong, burst of star for- 
mation may have occurred in the past, of which we now 
observe the intense extended radio component, and which 
could account for the overall steep non-thermal radio spec- 
trum of the galaxy. 

A realistic hypothesis is also that some SNRs and RSN 
are hidden because of strong free-free absorption, viz. in 
the 'mixed' regions mentioned in the previous subsection.^ 
Free-free absorption has in fact been also inv oked f or sev- 
eral compact sources in M82 (Wills et al. 1998), even 



though with turnovers at much lower frequencies, with the 
only exception of 44.01+59.6 that has a positive spectral 
index at frequencies higher than 1 GHz. 



^ For the unresolved sources synchrotron-self absorption that 
may lead to spectral indices up to 2.5 cannot a priori be ex- 
cluded. 



The combination of VLA and MERLIN 5 GHz observa- 
tions of the starburst galaxy NGC 2146 allowed us to re- 
solve seven of the eighteen hitherto unresolved sources 
detected with the VLA. We have measured the spec- 
tral indices between 1.6 and 5 GHz for nine sources. 
Three sources (37.6+24.2, 38.9+22.5, 41.4+15.0) seem to 
be good candidates for SNRs or RSN, with total emit- 
ted powers at 5 GHz comparable to the strongest SNRs 
found in M82. Four other sources (36.6+27.5, 39.0+17.9, 
40.2+18.1 and 37.7+23.7) have instead been identified 
with compact or ultra-dense H ii regions that have large 
EMs and a high-frequency turnover. The emitted thermal 
emission indicates a large equivalent number of 06 stars, 
exceeding the number of those found in similar objects in 
other galaxies. The large number of 06 stars may, how- 
ever, be consistent with stars in super starclusters. The na- 
ture of the last two sources (37.9+23.7 and 39.5+17.7) for 
which we have spectral index information are still a matter 
of discussion. It seems realistic that they emit a mixture 
of thermal and non-thermal radiation, and may manifest 
the most powerful H II regions ever observed. However, the 
burst of star formation in NGC 2146 has an extraordinary 
power, stronger than that in M82. Despite this and the 
fact that the overall radio spectral index of the galaxy re- 
veals a non-thermal origin, the number of discovered SNRs 
is much smaller than that in M82. For the explanation of 
this fact we propose two scenarios: 

— NGC 2146 is experiencing a burst of star formation 
stronger than that in M 82, but which we are observing 
in a younger phase. This phase still lacks a large num- 
ber of SNRs, but instead gives rise to strong thermal 
emission from optically thick dense ionized gas regions, 
associated with massive or super starclusters. The ori- 
gin of the non-thermal extended component must be 
a trace of past SN explosions, individually too old and 
faint to be seen, and belonging to a previous burst in 
the galaxy. 

— Strong free- free absorption by very dense Hii regions 
and, for the most compact sources, possibly also 
synchrotron-self absorption, may conceal a certain 
number of SNRs or RSN, rendering the identification 
of their nature difficult if only based on measurements 
of the spectral indices around the turnover frequencies. 
In this case the global number of SNRs or RSN can be 
higher, and this would corroborate the similarity be- 
tween NGC 2146 and M82. This interpretation raises, 
however, the question why we do not observe similar 
compact or ultra-dense Hii regions in M82 as well. 

Observations at 15 GHz will be able to better clarify the 
nature of the sources, and shed light on the unusual phase 
of the star burst we are observing in NGC 2146. 
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